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I/"") ■ A method of fundamental solutions has been used to investigate transitions in two 

energy level systems with no level crossing in a real time. Compact formulas for transition 
probabilities have been found in their exact form as well as in their adiabatic limit. No 
interference effects resulting from many level complex crossings as announced by Joye, 
Mileti and Phster (Phys. Rev. A44 4280 (1991)) have been detected in either case. It 
is argued that these results of this work are incorrect. However, some effects of Berry's 

^^ ■ phases are confirmed. 
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1. Introduction 

Transitions between energy levels in a two energy level system evolving in time are of 
great importance from many points of view. On one side such systems provide us with the 
simplest models to investigate transition amplitudes between different energy levels by differ- 
ent approaches Q. On the other side these systems play an important role in experimental 
investigations of basic principles of quantum mechanics Q. Recently a lot of effort has been 
devoted to obtain more rigorous results on the adiabatic limit of transition amplitudes for 
these systems || §, ||, [H fjj . In particular in a series of recent papers Joye et al have stud- 
ied this problem by the Hilbert space methods. Such two energy level systems are formally 
equivalent to a one-half spin system put into time dependent magnetic field. However good 
approximate results and the more so the exact ones are difficult to obtain for such systems 
even for simple time evolutions of the effective 'magnetic' field. Therefore each opportunity 
of improving this situation is worth trying. A treatment of the problem by a method of fun- 
damental solutions (so fruitful in its application to stationary problems of 1-dim Schrodinger 
equation [||, |j| |l(|) is of first importance, the more so that to our knowledge, the method 
was not used so far to this goal. A possibility of application of the method is related to 
the fact that a linear system of first order differential equations describing time evolution 
of transition amplitudes can always be transformed into a system of decoupled second order 
equations having a form of the stationary Schrodinger equation, one for each amplitude. This 



allows us to apply all advantages of the fundamental solution method [10]. The only obstacle 
related with this approach is a complexity of effective 'potentials' which appear in the final 
system of the Schrodinger-type equations. 

The paper is organized as follows. 

In the next section the problem of transitions in two energy level systems is stated and 
corresponding assumptions about the effective 'magnetic field' are formulated. A linear sys- 
tem of two differential equations for the transition amplitudes is rewritten in a form of two 
decoupled equations of the Schrodinger type. 

In Sec. 3 properties of the fundamental solution method are recalled. 

In Sec. 4 some subtleties of the application of the fundamental solution method to the 
problems considered in the paper are discussed. 

The method is first applied to a particular system of the atom - atom scattering within a 
frame of the Nikitin model |ll|, [l2| in Sec. 5 . 

In Sec. 6 results of Sec. 5 are next generalized to systems with an algebraic time dependence 
of the effective magnetic field. 

In Sec. 7 another two examples of two energy level systems are considered with corre- 
sponding magnetic fields depending exponentially on time. These examples show that a way 
the magnetic fields depend on time does not affect a form of the transition amplitudes. This 
form is not affected either by the number of (complex) energy level crossings on the Stokes 
lines closest to the real axis of the complex time plane. The latter result confirms the one 
of the previous section. Such a dependence resulting with some interference effects has been 
announced by Joye et al Q. 

In Sec. 8 we consider an example of the magnetic field with an explicit contribution of the 
geometrical (Berry) phase to the transition probability. 

We summarize and discuss our results in the last section . In particular we show there 
that the results of Joye, Mileti and Pfister [4] on the effects of interference from many level 
crossings are incorrect. 



2. Adiabatic transitions in two energy level systems 

In general, any two energy level system is formally equivalent to a one-half spin system 
put into an external magnetic field B(i). Therefore, we shall consider just such a system. 
Its Hamiltonian H{t) is given then by H(t) = ^fjJS(t) ■ a , where a = (<J x ,a y ,a z ) are Pauli's 
matrices so that two energy levels E±(t) of H(t) are given by E±(t) = ±^B(t) where B(t) = 

When the adiabatic transitions between the two energy levels E±{t) are considered then 
the following properties of the field B(t) are typically assumed to be: 

1° B(i) is real being defined for the real t, — oo < t < +00; 

2° B(i) can be continued analytically off the real values of t as a meromorphic function 
defined on some t-Riemann surface Re- A sheet of R_b from which B(i) is originally continued 
is called physical; 

3° On the physical sheet B(t) is analytic in an infinite strip S = {t : \Qt\ < 5,5 > 0}, 
without roots in the strip and achieves there finite limits for 3ft£ = ±00 , i.e. B(3ft£ = ±00) = 
B 1 * 1 7^ in the strip; 

The field B(i) depends additionally on a parameter T(> 0) i.e. B(i) = B(i,T) which 
introduces a "natural" scale of time to the system, so that its time evolution is expressed 
most naturally in units of T. If T is small in comparison with the actual period of the process 
considered then the latter is "fast" or "sudden". If, however, T is large in this comparison 
then the process is " slow" or " adiabatic" . 

In the adiabatic process of the system the following is assumed about B(£,T): 

4° A dependence of B(t, T) on T is such that a rescaled field B(sT, T) has the following 
asymptotic behavior for T — > +00 

B(sT, T) ~ BqOO + -B^) + ^B 2 (s) + ... (1) 

while its s-Riemann surface R^/T approaches 'smoothly' the topological structure of the 
Riemann surface corresponding to the first term Bo(s) of the expansion (|]). 

5° With respect to its dependence on s the field B (s) satisfies properties 1° — 3° above 
with substitutions t — ► s and B(s) — ► Bq(s). 

Note that condition 3° excludes periodic fields B(£) . 

The time-dependent Schrodinger equation induced by H (t) takes therefore a form 

^^f±= l -»B(sT,T)-ans,T) (2) 

The adiabatic regime of evolution of the wave function ^f(s, T) corresponds now to taking 
a limit T — > +00 in (|2|). 

The main problem of the adiabatic limit in the considered case is to find in this limit 
the transition amplitude between the two energy levels of the system for s — ► +00 under the 
assumptions that \l/(— 00, T) coincides with one of the two possible eigenstates ^±(—oo,T) 
of H(—oo) and that there is no level crossing for real t i.e. liminf B(t) > e > 0. Known 

— oo<£<+oo 

approximate solutions of this problem are that of Landau pl| and Zener p3] in a form of 
the so called Landau-Zener formula and that of Dykhne |l5| who have shown that such an 
amplitude should be exponentially small in the limit T — > +00. In the next sections we 
shall show how to get an exact (i.e. not approximate) result for this amplitude as well as its 
adiabatic limit with the help of the fundamental solutions . 



A typical way of proceeding when the adiabatic limit is investigated is using eigenvectors 
*f?±(s,T) of H(sT,T) satisfying (\P±, \I/±) = 0. Then, such eigenvectors ^£f±(s,T) can be 
chosen as the following ones 



^ + (s,T)=e- i /o^ sin2 f dCT 
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, ^_(s,T)=e- i /o* cos2 l d<J 



sin 
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where and <j> are polar and azimuthal angles of the vector B(i,T), respectively, and dots 
over different quantities mean derivatives with respect to s-variable. 
The wave function ty(s,T) can now be represented as 

*(s,T) = a+(s,T)e- iT ^ E +^ T ^qj + ( S: T) + a.(s,T)e- iT ^' E -^ T ^ d ^_( Sj T) (4) 

where s' takes any real but fixed value. 

The Schrodinger equation (g) can be rewritten in terms of the coefficients a±(s, T) as the 
following linear system of two equations 



a+ (s, T) = c(s, T)e { £ w(? ' T)d? a_ (s, T) 
.(s,T) = -c*(s,T)e- i /> (€ ' TK a + (s,T) 
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The system (|5|) can be rewritten further as the following linear system of second order 
equations 



(7) 



where the coefficient functions a± decouple from each other being however still related by 
By the following transformations 



+ ioj d-f + \c\ a + = 
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we bring the equations (||) to Schrodinger types 

b±(s,T) +T 2 q±(s,T)b±( S ,T) = 
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where 
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so that for real s (and T) we have 



Q-(s,T) = q* + (s,T) 



(11) 



The equations (|9|) are now basic for our further analysis since their form is just of the 
stationary 1-D Schrodinger equation. 

First let us note that the dependence of the "potential" function q + (s,T) on T is given 
by 



Q + (s,T) = \^ B i + ^ 
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where the dependence of -B, c, G, on T in (|i~2^) is also anticipated. By ( |i~2| ) we get a corre- 
sponding dependence of q-(s,T) on T. Taking into account ([j]) and (y) it is easy to check 
that the last formula provides us with the following type of asymptotic behavior of q+(s,T) 
for large T: 

q + (s,T) = $(s) + ^ (s) + ^ (s) + . . . (13) 

Therefore the above form of dependence of q±(s,T) on T permits us to apply to the 
considered case the method of fundamental solutions. For this reason we shall start the next 
section with a review of basic principles of the method suitably adapted to the considered 
case. 



3. Fundamental solutions and their properties 

Consider first q±(s, T) as functions of s. They are defined completely by an s-dependence 
of field B(Ts, T). According to our assumptions, the latter is meromorphic on some Riemann 
surface Rb/T. However, by (|12|), q±(s,T) are algebraic functions of B, B and B and, 
therefore, they are also meromorphic functions of s defined again on some other Riemann 
surfaces R-t determined by these algebraic dependencies. As it follows from ( |i"2| ) topological 
structures of H± can be quite complicated. However, in what follows, we are interested in 
the adiabatic limit T — > +oo by which the structure of R-t should be determined for T large 



enough basically by the first term q + (s) of the expansion (|T^). In consequence, by (|12|), it 
should be determined by /iB^°^(s) i.e. by the first term of the expansion (Q). The structure 
of R-t can turn out to be much simpler in this limit. Despite this supposed complexity 
of q±(s,T) and of their Riemann surfaces we shall introduce and discuss the fundamental 
solutions to the equations @ without simplifications. We shall do it for the q + (s,T) case of 
dl2|). An extension of the discussion to the q-(s,T) case will be obvious. 

A standard way of introducing the fundamental solutions is a construction of a Stokes 
graph ||, |9|, [l0| related to a given q + (s,T). Such a construction, according to Froman and 



Froman |8| and Fedoriuk ||], can be performed in the following way [10]. 

Let Z denote a set of all the points of R+ at which q + (s,T) has its single or double poles. 
Let 5(x) be a meromorphic function on R + , the unique singularities of which are double poles 
at the points collected by Z with coefficients at all the poles equal to 1/4 each. (In a case 
when R_|_ is simply a complex plain the latter function can be constructed in general with the 
help of the Mittag-Leffler theorem [17]. But for a case of branched R+ the general procedure 
is unknown to us). Consider now a function 

q + (s,T)=q + (s,T) + ±S(s) (14) 



The presence and the role of the 5-term in (14) are explained below. This term contributes 
to (|14l ) if and only when the corresponding 'potential' function q + (s,T) contains simple or 
second order poles. (Otherwise the corresponding 5-term is put to zero). It is called the 



L anger term [10, |lq| . 

The Stokes graph corresponding to the function q + (s,T) consists now of Stokes lines 
emerging from roots (turning points) of q+(s,T). Stokes lines satisfy one of the following 
equations: 



q+(Z,T)dt = (15) 

with Si being a root of q+(s,T). We shall assume further a generic situation when all the 
roots Si are simple. 

Stokes lines which are not closed end at these points of R+ (i.e. have the latter points 
as their boundaries) for which the action integral in (^) becomes infinite. Of course such 
points are singular for q+(s, T) and they can be its finite poles or its poles lying at an infinity. 

Each such a singularity Zk of q + (s,T) defines a domain called a sector. This is the 
connected domain of R + bounded by Stokes lines and Zk itself. The latter is also a boundary 
for the Stokes lines or being an isolated boundary point of the sector (as it is in the case of 
the second order pole). 

In each sector the LHS in (|15| ) is only positive or only negative. 

Consider now equation (0) for b + {s,T). Following Froman and Froman in each sector Sk 
having a singular point Zk at its boundary one can define a solution of the form: 

b+, k (s, T) = q-*(s, T)-e™ TW ^ T h + A^ T ) k = 1,2, . . . (16) 

where 

Ki r£.n-i 
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(17) 



with 



„ (S , TH ^L_W) + JL£?M) (18) 

ql(s,T) 4 ql(s,T) lb g*( Sj T) 



and 



W(s,T)= / Jq(£,T)d£ (19) 

where s, is a root of g(s, T) lying at the boundary of Sk- 



in fllq ) and (18) a sign of cr (=±1) and an integration path are chosen in such a way to 
have: 

aZ(W(&-W(t j+1 ))<0 (20) 

for any ordered pair of integration variables (with £o = s). Such an integration path is then 
called canonical. Of course, the condition ( p0| ) means that 6 +j fc(s,T) vanishes in its sector 



when s — > z^ along the canonical path. The Langer <5-term appearing in (14) and (|18| ) is 
necessary to ensure all the integrals in (fL8| ) to converge when Zk is a first or a second order 
pole of q(s,T) or when the solutions ([Lq) are to be continued to such poles. As it follows 
from ( |l8] ) each such pole z^ demands a contribution to 5(s) of the form (2(s — z^))~ , what 
has been already assumed in the corresponding construction of 5(s). 

4. The adiabatic limit in the fundamental solution approach 

Consider now the consequences of taking the large-T limit for the above description. 
We assume that for a given q + (s,T) and its Riemann surface R+ the corresponding Stokes 
graph G+ is drawn. It is drawn, of course, on the Riemann surface \/R+ corresponding to 

First let us notice that singular points of q+(s,T) such as its branch points and poles 
depend in general on T. For both kinds of these singularities this also means a dependence 
on T of jumps of q + (s,T) on its cuts as well as the T-dependence of coefficients of its poles. 

According to the property 4° of the magnetic field B (see Sec. 2) we can expect that 
the singular structure of q+(s,T), i.e. positions of its roots and poles, as well as the cut 
jumps and pole coefficients, change smoothly in this limit to their final positions and values 
respectively. This limit structure is defined by the singularity structure of q\ (s,T) (see 
expansion (|13|)). Therefore, both the topology of \/H+ an d the associated Stokes graph G + 

change accordingly to coincide eventually with the Riemann surface yR.+ and with the 

Stokes graph G^ corresponding to yq + (s,T) . This limit structure can be achieved in the 
following ways: 

a. some of branch points and poles of q+(s, T) escape to infinities of R+; 

b. some of branch points and poles of q+(s,T) approach the respective singularities of 

c. some of branch points and poles of q + (s,T) disappear because their respective jumps 
and coefficients vanish in the limit T — ► +oo. 

Being more specific we expect that for T large enough a set S+ of all singular points 
of q+(s,T) (i.e. containing all its branch points and poles) consists of three well separated 

6 



subsets S!l , S+ an and S+ . The points of S!T run to infinities of R+ when T — > +00. 
Those of S+ an disappear in this limit while those of S+ coincide in this limit with the set 
S^ of the singular points of q ( + '(s,T) . 

Let us remove the points contained in S+ U S+ an from the Riemann surface R + , i.e. 
let us consider these points as regular for q+(s,T). Then R+ will transform into R+ - a 
Riemann surface which singular points coincide with those of the set S^_ m . 

Together with the previous operation let us remove from \/~R+ a ls° the Stokes lines 
generated by the points of S!? U S'™ n so that the remaining Stokes lines can be uniquely 
continued to form the Stokes graph G+ generated by the set S^ m . It is clear that the graph 

G2 n coincides with G+ in the limit T — v +00. 

The above two operations will be called the adiabatic limit reduction or simply the re- 
duction operation. 

As we have mentioned earlier there is a set of sectors and a corresponding set of funda- 
mental solutions defined in them associated with the graph G+. By the reduction operation 
both sets can be reduced i.e. under this operation some sectors of G+ transform into cor- 
responding sectors of G+ whereas the others disappear. Obviously, the latter sectors are 
those which disappear when the limit T — > +00 is taken. 

The following assumption should stabilize the corresponding results obtained with the 
help of the fundamental solution method. 

6° Among a full set of fundamental solutions associated with the Stokes graph G + there 
is a subset of them associated with graph G+ which allows us to solve the basic problem of 
the adiabatic transition and which is invariant under the reduction operation. 

The dynamical systems described by the Hamiltonian H(t) satisfying assumption 6° will 
be called the adiabatic limit reducible (ALR-)systems. 

The above assumption means that to solve the problem of the adiabatic transitions in the 
ALR-system we can first perform the reduction operation and next work with the simplified 
Stokes graphs G+. A set of fundamental solutions associated with this graph which can be 
used to solve the problem considered coincide with the corresponding ones of the full graph 
G_|_. The procedure used to construct a solution of the problem with the help of the latter 
graph is not affected by the reduction operation, i.e. it looks the same when the simplified 
graph G+ is used instead of G+. Therefore the aim of the reduction operation is to make 
easier choosing the proper set of fundamental solution solving the problem. The results 
obtained in this way can be still exact if the integration paths taken on the graph G+ can 
be mapped properly on the Stokes graph G+ restoring in this way the exact condition of the 
problem. However, if such a map is not known or is difficult to construct (because of the 
complicated structure of graph G+) the result obtained in this way can be considered only 
as an approximation i.e. valid only in the limit T — * +00. 

According to the above assumptions we can conclude from ( Jl2[ ) and fll3| ) that there is 
one-to-one correspondence between the Stokes graphs G + and G+ and the corresponding 
sets S+ and S+ . Namely, this correspondence is built by aggregations (blobs) of singular 
points of S+ , i.e. the branch points and poles of q+(s, T), which are transformed into single 
points of S + when the limit T — ► +00 is taken. Also there are sheaves of Stokes lines of 

fin (0) 

Gi. emerging from the blobs and transformed into single lines of G\ in the same limit. 



Therefore in the limit T — > +00 we can eventually consider for potentials (|12| ) Stokes 
graphs corresponding to first terms q± (s) of the asymptotic expansions for q±(s,T). The 



first terms of the asymptotic expansions corresponding to q± (s) and q±(s,T) are the same 
in this limit and equal, according to ([!]), to ^ / u 2 Bo(s). 

Let us note that properties 1° — 6° above can be satisfied by the field B for which B 2 
is a meromorphic function of t. We shall assume just such a dependence of B on t and of 
the corresponding rescaled field B(sT, T) on s. However, for simplicity, instead of continuing 
our considerations in their most general form we shall investigate first a particular example 
of the field B(t, T) which, as it seems to us, will illustrate our method in a satisfactory way. 



5. The Nikitin model of the atom-atom scattering 

The model of Nikitin [O] describes the scattering A*+B— >A+B+Ae of the exited atom 

A* moving with a small velocity v with the impact parameter b' and scattered by the atom 

B. The interaction between the atoms is of the dipol-dipol type. The latter example was 

analyzed in the context of the adiabatic limit v — ► also by Joye et al Q . 

The Hamiltonian for this system reads (PH, paragraph 9.3.2 and p 
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where Ae and C are constants and R = v b' + v 2 t 2 is the distance between the atoms. 
Introducing d = (2C/Ae)3 as a natural distant unit for this case and T = d/v as the 
corresponding adiabatic parameter and rescaling: t — ► sT and b' — > bd we get from (pl|): 
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In the 'magnetic field' language we have of course B(sT, T) = I (b + s 



3 

2 ,0,1 



so 



Ac 
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problem is determined by the function 

The graph is shown on Fig.l. 

Each g-t (s, T) has 40 roots, five branch points at s = ±ib and at s = s& = ± I e 3 ~b ) 

, A; = 1,2,3, as well as two poles at s = 0. Therefore only six roots of q(°'(s,T) at s = Sk, 
A; = 1,2,3 and its two poles at s = ±ib look encouraging. Nevertheless, we shall consider 
first the case without any approximations. 

At first glance the Stokes graphs corresponding to the functions q±(s, T) seem to be quite 
complicated. However it can be handled in the following way. 

Functions q±(s, T) are determined on two sheeted Riemann surfaces R-t respectively with 
the branch points at s = ±ib and at s = s^, k = 1,2,3 and with 40 roots distributed into 
halves on each sheet of the surfaces. Therefore the Riemann surfaces vR± corresponding to 
y/q±(s,T) (it will turn out that it is not necessary to introduce to the latter functions the 
corresponding Langer terms) are four-sheeted with these 40 roots being square root branch 
points on them. When T — ► +cx) only six of these branch points survive coinciding with the 
six roots of q(°'(s,T) at s = ±s/«, k=l,2,3 whereas H± transforms into the complex s-plane 
since the branch points of q±(s,T) at s = ±ib disappear, being transformed into the second 
order poles of q(°'(s,T). It is easy to check however that for finite but large T these six 
roots of q(°'(s,T) are each split initially into two as. The split is the result of the square 
root branch points at s = ±ib to which the recovering of the finite T transforms the poles 
of q(°>(s,T) at the same points. The two copies of each of these six roots lie of course on 
different sheets of R±. Next, each of these 12 roots is still split into three by the same reason 
of finiteness of T. In this way, on each of the two sheets of R-t there are 36 roots grouped 
by three around their limit s = ±s/c, k = 1,2,3 achieved for T — ► +00. 

The remaining four roots of q±(s, T) are displaced in two pairs, one pair on each sheet of 
R-t, close to the points s = at which the second order poles of q±(s, T) are localized. When 
T — ► +00 the roots in each pair collapse into s = multiplying the corresponding second 
order poles and thus causing mutual cancellations of the latter and themselves in this limit. 

Now we shall focus our attention on the Stokes graph G_ generated by q-(s,T) on 
the first sheet of R_ as well as on the remaining ones. It looks as in Fig. 2. The Stokes 
graph G + corresponding to q+(s,T) can be obtained from G_ by complex conjugation of 
the latter. On the figure the wavy lines denote the cuts corresponding to the branch points 
of the fundamental solutions defined on R_. The sheet on Fig. 2 cut along the wavy lines 
defines a domain where all the fundamental solutions 6_ : i(s,T), ...,&_ 2(5, T) defined in the 
corresponding sectors Si, ..., S^ (shown in the figure) are holomorphic. 
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Fig.l The Stokes graph corresponding 
to 'potential' (M) 



Fig. 2 The Stokes graph corresponding 
to 'potential' q-(s,T) of (23) 



According to our earlier description of the behavior of the Riemann surfacey r R+ when 
T — ► +00 the set S_ corresponding to the considered case is empty, S™ n contains four 
points at s = on each of the four sheets of ^/R_ (these four points correspond to the second 
order poles of q-(s,T)) and the four branch points close to s = 0, while S„ contains all the 
remaining singular points of \/g_(s,T). Now, for our case, the solution of the problem stated 
in this paper is simple. Namely, it can be found in the following steps: 

i. take a linear combination of the fundamental solutions 6_ : i(s,T) and b_i(s,T) to 
construct the amplitude a_(s,T) with the desired property at s = —00, i.e. lim |a__(s,T)| = 

s— »— 00 

0. This amplitude is defined in this way up to a multiplicative constant; 

ii. use the equation (||) to construct a + (s,T) and adjust the constant mentioned earlier 
so that the limit lim |a+(s,T)| = 1 can be satisfied; 

s^— 00 

iii. continue canonically a^(s,T) along the real s-axis with the help of the solutions 
b- t i(s,T) and b_ i(s,T) using to this goal the remaining fundamental solutions if necessary; 
iv. calculate the limit s — ► +00; and 
v. calculate the adiabatic limit T — > +00. 
According to @ and to the first of the above steps we have: 

a-[s,T) = Aq_ i (s,T)e Js v Js o xi(s,Tj 



(25) 



+Bq_ 4 (s,T)e Js ' 2Kc Jy ' Js H -^ ' J Xi(s,T) 



where s' is any point on the real axis which is regular for the integrand whilst sq is the one 
from the infinite strip bounded by the Stokes line abode from one side and by a'b'd from 
the other (see Fig. 2), being also an arbitrary but regular point for all the integrands. The 

choice of signatures in (p6|) was done due to the fact that 3? I iT J s s qlda J is positive (for s 

sufficiently large) for the sector Si and negative for 5j . The latter property follows from 
the fact that according to d2^) and the Stokes graph on Fig. 2 we have on the first sheet of 
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^/R_: sgn Ky ql(s,T) =sgn(s) for s — ► ±00 along the real axis. 

If further we take into account the following asymptotic behavior of the relevant quantities 
on the real axis: 



1 


c 




2 


c 


- IUJ 


1 


"c 


1 


2 


_c 


- ZCJ 



+ iTJq- 



-iTAe 



iT^T~0(- 



-00 



-00 



(26) 



then we can conclude that B = in (Pq). 

To fix the value of the constant A in (^) we can use the second of relations (|5|) and 

apply the condition mentioned in the second step of the procedure i.e. lim 



i I , uida ■ 
' J s' ( 



a-(s,T) = 1 to get: 



A 




exp 



*o 



iujds + 



*'<) 



TAeV 2 
Therefore, for the amplitude a~(s,T) we obtain finally: 



iuj + iTy/qi 



ds + In c(sq] 



c(s,T) 



(27) 



a-(s,T) 



1 /Ae -i 
TAlVl"' 7 - 4(S ' T)eXP 



«o 



iiods 



"2 (c~ iu; ' ! + ' :/ v ^ 



ds+ 



lnc(so) + / 

■/sn 



(28) 



i(£-ia;]+,:/\ 



rfcr XI (s,? 1 ) 



Now we can take the limit s — ► +00 in the above formula, continuing along the canonical 
path 7i^ 2 shown in Fig. 2, to get: 



a_(+oo,T) 



1 



iTAe 



exp 



lnc(so) + 



«o 



*o /-so 

icjcis + / 



1/c . 
2 Vc 



1/c . 

2 Vc 



+ iT, 



ds+ 



(29) 



+ iTy/ql 



d^ XI-^ 2 ( T ) 



The apparent so-dependence in the above formula is illusive. We can use this fact to 
calculate the integrals in the exponent most accurately. First let us note that we cannot 
disjoint totally the integrations in the two infinite integrals since the diverging contributions 
of the three terms in both of these integrals cancel mutually at the corresponding infinities, 
making the integrals convergent. We can however take as the integration paths for these two 
integrals the Stokes lines abc on Fig. 1 and abode on Fig. 2. Namely, let the points sl on 
line a and sr on line e be arbitrarily close to the corresponding infinities of the real axis. Let 
further points s' L and s' R be the points on the Stokes lines a and c of Fig. 1 respectively. We 
choose the latter points to lie on the anti Stokes lines of Fig. 1 which pass by the respective 
points sl and sr. Then the integral in the exponential of formula (|30|) can be rewritten as: 



+ iTJq~: 



SO 


/■•so 


r 1 


fr 




iojds + 


/ 


-- 


[-- 


- ILO 


1 j 


'-00 


2 


\c 





ds + lnc(so) 
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+ 



+ 



1 ft 



+00 



*() 



i(!-*"i + ,T - 



ds 



iujds 



2 Vc 



+ 



+00 



1 /c 



2 Vc 



iw + iTJqi 



iuj ) + i7\ 



1 f S L 1 

ds + - iuds + - In c(sl) 

2 Js L 2 



1 /" s k 1 

+ - / iujds + - In c(s R ) 

4 Jsn ^ 



(30) 



rs R 1 

+ / iT^fqZds + - 

Js L 2 



iuods 



iuods 



where the last three integrals run along the respective Stokes lines and, therefore, are purely 
imaginary. Point s" in the above formula is an arbitrary point of the Stokes line 06c on Fig. 
1. 

We are interested mainly in the transition probability defined by amplitude a_(+oo,T) 



for which only the real part of the integral I is important. Formula (30) gives for it: 

1 [ s 'l 1 

ds + — I iuods + -3? In c(sl) 

2 J St. 2 



til- -T« /" icuds + M I ' 

J — 00 



+3? 



■*•/? 



"2 (~ ~ iuJ ' + ll v ' /_ 



1 /c 



(31) 



2 Vc 



icj + iTJTf- 



1 /" s fl 1 

ds + - I iuods + —3? In c(sh) 



*i? 



We can now calculate 3?I taking in (|32|) the limits s^ — ► —00 and s/j — > +00 along the 
corresponding Stokes lines. We get in this way: 



3tt = -3? / iwds + - lim (/ L iuods + 3? In c(sx,) ] + 
Js' 2 s L ^-oo yy SL y 

- lim / iuds + 3? In c(sr)\ =—3?/ iuods + In - 
2 s R -*+oo yj SR J j s , 2 



(32) 



The limits in ( |32[ ) can be obtained by estimating the asymptotic behaviour of the differ- 
ences s' L R — sl,r and the corresponding functions when \s\ — > 00 along the Stokes lines, for 
which direct calculation gives: 



s l,r ~ s l,r 



4i 

TAe 



In I s I 



5b 2 i ihia^R 



2TAe\s\ 2 TAe 
iuj ~ iTAe ( 1 + 



(33) 



2s 6 

2 56 2 

3? In c(s) ~ In 4 In I si t—tk 

K ' 3 " 2|s| 2 

where constants o,l,r are also independent of T and can be estimated exactly only when the 
exact equations of the Stokes lines 06c of Fig.l and abode of Fig. 2 are known. 
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The imaginary part of the integral I can be calculated as the following limit: 

7(T) = 9fI= lim 9 ( -In c(sl) + -In c(sr) 
slm^T^o \2 2 

1-sr ^ rs" \ rs' R rs" \ 

+ / iT^/qlds H — / iu>ds / iuods — / iuods (34) 

Js L 2 7 S ' L 2 7 S " 7s' / 

Therefore, the final exact formula for the transition amplitude is: 

I i rp\ 3 a LO-R - f s ; iui(s,T)ds+i-y(T) t rp\ / Q r\ 

a_(+oo,T) = 2TAg e Js ' ; Xi^ 2 ( T ) (35) 

and the probability P reads: 

4T2(Ae)2 e IXl-» 2 (r)| (36) 

where in the last two formulae point s' is an arbitrary point on the real axis while point 
s" being the one of line abc of Fig. 1 is taken to lie simultaneously on the anti-Stokes line 
passing by point s'. 

The adiabatic limit of the transition probability is therefore: 



_ 9aW R „- m ( iT f*> » B °(s)ds) 

^"4T2(Ae)2 e [67) 

where s" is now an arbitrary point of the continuous Stokes line passing by roots of Bq(s) 
closest to the real axis. 

6. The general case of algebraic magnetic field 



The result given by the formula ( |35|) can be easily generalized. From the way of obtaining 
formula ( |30|) it follows that the most important is the existence of the continuous Stokes line 
abode on Fig. 2 and its T — ► +oo-limit, i.e. the Stokes line abc of Fig. 1, which link the 
respective infinities 3fts = -co and 5R = +oo on both Stokes graphs. Another important 
property was the way field B approached the limits B when 3ft£ — ► ±oo respectively in the 
strip 5] mentioned in the assumption 3°. Let us therefore accept the following two additional 
assumptions: 

7° There are two Stokes lines on each of the Stokes graphs corresponding to iT^Jq± 
which can be taken as the boundaries of the strip X. Each of these two Stokes lines links 
continuously both infinities of the strip S, see Fig. 3; 

8° Inside the strip £ the field B approaches the infinities of the strip according to the 
following asymptotic formula: 



B( S T,T)~B±(T) + 5£p + 5iip + ... + 5£p + ... ) Ks ^ ±00 



S a l S a 2 s"fc 



1 ± ± ± 

- < Of < aj < • • • < a k < ■ ■ ■ 



(38) 
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where a 1 , . . . , a k , are rational if B 2 is a meromorphic function of s. 

If the Stokes graph corresponding to iT^fqZ satisfies the conditions of being a graph of the 
ALR-system described in Sec. 4, then we can claim that there are four sectors Si, £j, £2, S^, of 
the graph and the corresponding fundamental solutions xi>Xi which can be used in exactly 
the same way as it was done in the case of the Nikitin model to solve the problem stated in 
Sec. 2, see Fig. 4. 




Fig. 3 The Stokes graph corresponding 
to general q(°'(s) considered in Sec. 6 



Fig. 4 The Stokes graph corresponding 
to general q-(s,T) considered in Sec. 6 



Let us choose the xyz-axes in the space of vector B in such a way that one of its limit 
components B x and B does not vanish in the corresponding infinities. Let us also assume 
that vectors B (T) and B x (T) of expansion (^8|) are not parallel to each other in the respec- 
tive infinities (otherwise we should take another pair of vectors appearing in (^) satisfying 
the last property and having the smallest sum of the power exponents by which they are 
accompanied). Then, if we take into account the following asymptotic which comes out of 
(ph and of the above assumptions when Ks — ► ±00 inside the strip: 



B±x 



B, 



±2 



uJ 



^TB± 



B£ x B± 



r?=L2 1 r>±2 

B„ n + B yQ 



'x,0 



I 

+ 2 



BjxB± 



ni / E?i2 1 r>i2 



M TB± • B± 



B 



z,a 



(B± x Bf 



m 



rj=b2 I ni2 
n x,0 "1* n y,0. 



s a i 



s a t 



D± 



s a i 






1 
2 

1 

2 



iuj 



IUJ 

c* 



+ iTJq~- 



iTyfqZ 



-ifiTB^ 

p-(D-_r 1 

D+jD+T 1 



(39) 



,± 



o 1 
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then we can repeat the procedure of the previous section to get the analogues of formulas 
( |30| ) and (gj). Namely, we have for them: 



a_(+oo,T) 



1 



VTJB {T)B+(T) 



■. cxp 



«o 



lnc(s ) 



so 



iujds + 



2 \c~ iU) ' ] + '' / ' v '' / " 



w ) + iTJqi 



ds 



d<? X1-+2OO 



(40) 



ol^r / |£>-(T) J p+(T)i c _ r ; MSiT)ds+ , 7 

/iT V Bq(T)B+(T) € 



XI^2(T) 



where points s' and s" have been chosen again on the same anti-Stokes line of the graph 
corresponding to iu>(s,T) and 



_ alal\D-(T)D+(T)\ mr ;' MsiT)d 2 

F{T) - (^TYB Q {T)B+(T) e "^ IXi-2(r)| 



(41) 



where D^ are given by: 



so that: 



D d 



B± x (B± x B± 



R ±2 



r?±2 I r?±2 
B x,Q + ^y,0 



+ 



(B? x B±), 



Bn\/B±Z + £ 



,±2 



(42) 



l-D 



±| __ S^sin^ 
2B£ 



(43) 



where ^(T) are the angles between fields B and B x respectively. 

Again, the exact form of the coefficients cll,r can be found if the exact equations of the 
Stokes lines corresponding to u>(s,T) and g_(s,T) are known. 

Therefore, the final forms of the transition probability and its adiabatic limit are: 



pm _ g|4gj~ (T) Bj (T) sin 0- (T) sin 0+ (T) 2SR r" Ms ,T)«fo, m|2 

(2^TB^(T)B+(T)) 



(44) 



and: 



jad 



a z L aj i B 10 B{ sm(j) sin^J -2K(iT/ s s , )j,B (s)ds 



(2 M TS - S + 



(45) 



where to get the last formula the asymptotic expansion (ffl) has been applied to fields B (T) 
and B x (T) as well as to ui given by (||). Point a" is now an arbitrary point of the continuous 
Stokes line a\a2 ■ ■ ■ a n a n+ \ passing by the roots of Bq(s) closest to the real axis, as it is shown 
on Fig. 3. Note that because of our assumption the angles in ( f44[) and (^5|) are different from 
and tv. 
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7. Another two examples with exponentially decreasing mag- 
netic fields 

We consider here another two examples of magnetic fields depending exponentially on 
time. The main difference between these cases and those considered in the previous sections 
lies in the number of level crossings which in the exponential cases is, of course, infinite. 
We consider the following two cases of the fields: 

a) B(t,T) = B + — -^y, B • Bi = 0, B = \B \ ^ \B X \ = B ± 

cosh [7p) 



(46) 



b| B(/.7-)-B„ + B; L lanli(|;). B n -Bj-0. |B„| - B, | - B {] 



Case a) 



The relevant quantities for this case have the forms: 

_ 1 BpBisinhs 



2B( + B$ cosh 2 s 



u{s,T)= f j,T\Bl + ^- 
V cosn s 



2 



, , BZsinh(2s) _ / , B? , 

( S T) = -— - coths =-5 — - v ; 9 -inT\ BZ + \-\ + 

4T 2 I B\ + £ 2 cosh 2 s P V cosh 2 s ' 



(47) 



2,.;«l,')„ 1 r2d2 c ,:„i 1 2 



i/x 5f sinh 2s 1 ^Sf sinh 2 s 



r + 



2T • ^,2 B ( R 2 , R 2 _u2 \ 2 4T 2 / „2 , R 2 ms u2 _ 



sinh 2 s 5 2 + £ 2 cosh 2 s 2 ^ LB 2 + 5 2 coshes 



2B 2 cosh 2s B# sinh 2 2s 



2T 2 



Sinh 2 S 5 2 + Bl COSh 2 S fn2 , n2 u2 n ' 



ysmii s £>! f ^ cusn s ^ B z + B g cogh z s 



and the Stokes graphs defined by w(s, T) and q-(s, T) are shown on figures 5. and 6. respec- 
tively. 
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|S T 



fs 2 



Res 




\S- 2 




Fig. 5 The Stokes graph corresponding 
to q(°' (s) of case a) of 



Fig. 6 The Stokes graph corresponding 
to q~(s,T) of 



The procedure leading us to formula ( |30| ) is still valid but the corresponding sectors 
Si, Si, 5*2, 5*3 are now less exposed. Namely, the first two lie on the left of the imaginary 
axis, Si above and Si below the real axis whereas the next two on lie the right of the imaginary 
axis and, respectively, above and below the real axis. A peculiarity of this and the next case 
is that these sectors are cut by the infinite number of the Stokes lines parallel to the real 
axis and distributed up and down to the imaginary infinities, see figures 3. and 4.. The 
fundamental solutions defined in these sectors vanish in their imaginary infinities. Therefore, 
the corresponding transition amplitude a,-(s,T) from level E+ to E- looks as follows: 



a-(T) 



2iBiaLCtR 



HB Q JI + ^T^BI 



■jfMs,T)ds +l y(T) x _^ T ^ 



(48) 



To get the above formula we have taken into account the following asymptotic behaviour 
of the quantities determining it: 



c(x + iy) 



+ 



B 



i p -x-vy 



Bo' 

Bxp+x+iy _ 
B C X 



+00 

— oo 



u(x + iy, T) ~ -fiTB 



oo, 



Vl,r ~ Vl,r 



c(s) r 

2\na L)R \xl,r 



-1 Ks -> +oo 

+1 JRs -^ -oo 



(49) 



+ 



fjTB fJ<TB ' 



\ x l,r\ 



X 



where sl,r = xl.r + Wl,r and s' L R = x' L R + iy' L R have the same meaning as previously, 
i.e. lie on the corresponding Stokes lines defined by q_(s,T) and lj(s,T), respectively, whilst 
cll,r measure (together with the terms linear in xl,r) the deviations of these lines at the 
corresponding infinities. 

Therefore, for the exact transition probability and its adiabatic limit, we obtain from 



P{T) 



(2Bia L a R ) 



tPBKl + n^Blf 



-2R/;"«*(«.:n<fo| XI _ +2 ( ;r )| l 



(50) 
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and 



respectively. 



Case b) 



P 



- { fTBl ) 6XP 



-2/jTU 



«/5g + 






cosh s 



(51) 



In this case we have: 
c(s,T) = - 



1 



2cosh(2s) 



c(g,r) 

c(s,T) 



-2tanh(2s) 



w(a,T) = nTB 



\/cosh(2s) 
coshs 



(52) 



«-(s,T) 



1 / , / n ^^ \/cosh(2s 

' 2tanh(2s) + ifiTB 



4T 2 
i/xBq tanh(2s) — tanh s 



2T 



cosh s 



cosh(2s) 



cosh s I 
7 1 



4T 2 cosh 2 (2s) 

and the Stokes graphs corresponding to u>(s,T) and q^(s,T) are shown in figures 7. and 8. 
respectively. 



i, lm s 



Res 




Fig. 7 The Stokes graph corresponding Fig. 8 The Stokes graph corresponding 

to g(°) (s) of case b) of © to g_ (s, T) of 



Again, the transition amplitude can be calculated taking into account the following asymp- 
totic: 

,-2x-2iy x _, +00 ^ ( _ 2 Sfo _+ +QO 

-oo c(s) 1 +2; Us — >■ -oo 



c(x + zy) 



_ e +2a;+2iy ^ 
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uJ 



(x + iy, T) ~ — y2fiTBo, \x\ — ► oo, (53) 



, 2lna LtR xx, h 



so that we get for it: 



_( T ) = 2 ^ Qfi e ~/; «(.,T)«fe-H7(T) (r) (54) 



Therefore, for the corresponding transition probabilities we obtain: 

D/rpN, (2axO/j) _ 2 R r i w ( s T)ds| //riM 2 / KK \ 

P(r) = ^(4 + 2^r^g) e Js lx ^ 2(T)l (55) 



and 



pad (\/2a L a R \ I f s " . /cosh(2s) ' 

p = b^j exp \-^ TB ^j sl i^^ d \ 



(56) 



8. Non vanishing contribution of the Berry phase 

The previous sections have provided us with the examples of Hamiltonians in which the 
corresponding transition probabilities have had no contributions from the term 

of to (see (0)) representing (at least) a part of the Berry phase of the transition amplitudes. 
The Hamiltonian defined by the following field: 



B z (BxB 
' B Bj+B$ 



B 



B n 



1 + s 2 



1, as, s 



a > s/2 (57) 



provides us with the corresponding positive example of such a contribution. Namely, for this 
case we have: 



v (l + s 2 ) 2 + a 2 -2 l s 2 

uj = fiTB ——2 — — 5-5 (58) 

1 + 5 /(i + s 2f + Q 2_ 2 l + aV 



From (|57| ) it follows easily that for this case the transition probability (|44| ) takes the form 

(2^T£ ; 



P(D = /n a ^R 2 e" 2 *//' " (s,T)ds |xi_2 (^)l 2 (59) 



It is the second term of (|5g) which is responsible for the Berry phase contribution to 
the transition probability (|59|). We shall calculate this contribution in the adiabatic limit 
only and for a close to \/2. This assumption allows us to calculate the corresponding path 
integral: 



s 2 



Iv = — i / . =- — „ ds (60) 

A J(l + a 2)2 + Q 2_ 2 l + a2 S 2 
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^ 2 + a 2 - 



from point s = to the closest root sq = iyl + iyor — 2 of the polynomial (1 + s z y + a l — 2, 
lying in the second quadrant of the s-plane. For a close to y2 we can simplify the integration 
expanding suitably the square root in the integrand of ( p0| ) and the root sq as well. It is easy 
to check that under the above assumptions the net result of such calculations is: 



23HL = In ■ 



a 



V2 



22 



*(72-l) 



v/2 



+ o 



a 



V2 



(61) 



Obviously, the above Berry phase contribution to the transition probability ( p9[) modifies 
its prexponential factor multiplying it by the following additional one: 



C 



a 



V2 



22 



2(^2- 



v/2 



(62) 



9. Conclusions and discussion 

We have shown in this paper that the fundamental solution method has turned out to 
be very effective also in its application to the problems of the transition amplitudes in two 
energy level systems. In particular it has enabled us to obtain compact and exact formulae 
for these amplitudes and to get easily their adiabatic approximations as well. Due the clear 
way of their obtaining and their compact forms, the formulae allow us to claim that there 
are no particular effects coming out of the many complex level crossings i.e. there are no 
of individual contributions of any kind to the transition amplitude from each such crossing 
leading to any particular interference effects in these amplitudes. Just the opposite, such a 
contribution is controlled totally by the Stokes line closest to the real axis of the t-plane which 
is however built by these crossing points of the two energy levels. This result is independent 
of both the number of complex level crossings (i.e. finite or infinite) and of the particular 
type of the t-dependence of the effective magnetic field (i.e. algebraic or exponential). In 
this way the respective results of Joye, Mileti and Pfister [4] have not been confirmed by 
our approach. This last difference seems to be rather dramatic and, as it seems to us, its 
origin lies in an erroneous calculation of the transition amplitude by the authors mentioned. 
Namely, it is the formula (6.21) of their paper [4] for the transition matrix X{z\) which is 
wrong, particularly if applied further in their Lemma 6.1 to get the general formula of it. 
This can be seen if we rewrite our results in terms of the transition matrix. 

Namely, let us denote the r.h.s. of the formula (28) by U2i(s, T) and the result we obtain 
calculating a + (s,T) with the help of the second of Eqs.(5) and of (28) by Un(s,T) . 

Reversing the problem we have solved in our paper by assuming that for s = — oo the 
vanishing amplitude is rather a+(s, T) than a_ (s, T) we obtain by exactly the same methods 
as used in Sec. 3 and the further ones the remaining elements Ui2(s,T) and U22(s,T) which 
construct the transition matrix U(s,T). For the choice we have done in our paper we have 
of course: 



a + (s,T) 
a-(s,T) 



U u (s,T) U 12 (s,T) 
U 2 i(s,T) U 22 (s,T) 



a(«,T) = 

Matrix U(s,T) is of course unitary (for real s) and U(— oo,T) = I 



(63) 
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It should be clear that the order of Ui2(s,T) as a function of its arguments is the same 
as that of U2i(s,T), the latter element being given by the r.h.s. of (28) so that the adiabatic 
limit of U(s, T) is given by: 



U ad (s,T) 



Ug(a,T) 



U$?(s,T) 1 



(64) 



Therefore this matrix is not a triangular one in this limit, as it is the case of X(z\) 
mentioned earlier, which does not contain the non vanishing element U^{s,T). 

Moreover we can not apply matrix U(s, T) directly to continue the solution (p3|) along the 
central strip of the corresponding Stokes graph from point s to another one s' . The proper 
continuation is of course the following: 



a( S / ,T) = U( S , ,T)U- 1 ( S ,T)a( S ,r) 



(65) 



In particular, if it is possible to continue the solution a(s, T) along, say, the upper Stokes 
line limiting the central strip (i.e the level crossings si, S2> ••• j s n met along this line are not 
an obstacle to such a continuation) then continuing a(s,T) in this way to s = +oo we get: 



a(+oo, T) = U(+oo, T)U- X (s n , T)XJ(s n , TjlT 1 (s n _i , T) .. 



(66) 



\J(s 2 ,T)V- 1 (s 1 ,T)\J(s 1 ,T) 



U(+oo,T) 



The above results show that none of the contributions from the individual level crossings 
lying on the considered Stokes line survive on the way of continuation. 

On the other hand writing both formula (6.21) and the respective result of Lemma 6.1 of 
[4] in terms of the quantities introduced above, we get: 



*o»i: 



i 







Ug( 8l ,T) 1 



(67) 



and 



a ad (+oo, T) = X{s n )X{s n ^) . . . X(s 2 )X( Sl 



1 




Comparing the last two formulas with the respective (|64]) and (|6^) ones we see that 
formulas ( p7| ) and ( p8| ) are wrong. Particularly, it is the incorrect formula (B3) which gives 
rise to the interference effects in the amplitudes of Joye, Mileti and Pfister [4]. 

Finally, we would like to mention that, as we have shown this in Sec. 8 there are con- 
tributions to the transition probabilities originating from the geometrical (Berry) phase [16] 
although their geometrical meaning in the context of the transition amplitudes is not clear. 
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